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Abstract 


Thermodynamic analysis of dimethyl ether steam reforming (DME SR) was investigated for carbon formation boundary, DME conversion, 
and hydrogen yield for fuel cell application. The equilibrium calculation employing Gibbs free minimization was performed to figure out the 
required steam-to-carbon ratio (S/C =0-5) and reforming temperature (25-1000 °C) where coke formation was thermodynamically unfavorable. 
S/C, reforming temperature and product species strongly contributed to the coke formation and product composition. When chemical species DME, 
methanol, CO2, CO, H2, H2O and coke were considered, complete conversion of DME and hydrogen yield above 78% without coke formation 
were achieved at the normal operating temperatures of molten carbonate fuel cell (600°C) and solid oxide fuel cell (900°C), when S/C was at or 
above 2.5. When CH, was favorable, production of coke and that of hydrogen were significantly suppressed. 


© 2006 Elsevier B.V. All rights reserved. 


Keywords: Thermodynamic analysis; Carbon formation boundary; Dimethyl ether; Steam reforming; Fuel cell; Hydrogen 


1. Introduction 


Fuel cells are presented as an efficient power generator appli- 
cable to both mobile and stationary use. In typical, an electro- 
chemical reaction of a fuel at anode (generally hydrogen) and 
an oxidant at cathode (generally oxygen) generates electricity 
through a fuel cell. Several types of fuel cell have been devel- 
oped, for example, polymer electrolyte fuel cell (PEFC), solid 
oxide fuel cell (SOFC), and molten carbonate fuel cell (MCFC). 
PEFC provides high power density and operates at low temper- 
ature of ca. 80°C. These features make it suitable for vehicle 
application. However, carbon monoxide in the fuel feed is lim- 
ited generally below 10 ppm since it poisons Pt anode of the fuel 
cell, especially at such low temperature. Therefore, pure hydro- 
gen or external fuel reformer with CO remover is needed for 
PEFC. Unlike PEFC, SOFC and MCFC operate at high temper- 
atures of 700—1000 and 600-650 °C, respectively. Both external 
and internal fuel reformers are applicable to these types of fuel 
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cells. In the case of the internal reforming, hydrogen is produced 
by the reforming reaction taking place on the fuel cell anode 
or inside the fuel cell chamber. In addition, intermediate tem- 
perature fuel cells operating at 400-600 °C has been presently 
developed. 

Various types of hydrogen production systems have been 
extensively studied. Hydrogen is an environmentally friendly 
fuel since it does not release carbon dioxide on site when 
employed in fuel cell processor. Several hydrogen-generating 
techniques, i.e. steam reforming, partial oxidation, and autother- 
mal reforming of various fuels, e.g. gasoline, LPG, methane, 
methanol (MeOH), ethanol (EtOH), and dimethy] ether (DME), 
have been regarded as the efficient processes for deployment of 
fuel cell systems. Presently, DME is recognized as a promising 
alternative hydrogen source. Of DME synthesis, MeOH is first 
synthesized from syngas over metal catalysts for example Cu- 
based catalysts [1-4]. Subsequently, dehydration of MeOH to 
DME takes place over solid-acid catalysts for example zeolite 
and alumina [1—4]. Several large DME plants are currently under 
construction [5]. The large scale production will decrease the 
cost of DME even though the cost of DME is still controversial. 
DME and MeOH are suitable for on-board reforming because 
they can easily be liquefied and be reformed at low temperature 
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Nomenclature 

dik number of atoms of the kth element present in 
each molecule of species i 

Ak total mass of kth element in the feed 

fi the fugacity of species i in system 

1 the standard-state fugacity of species i 

FDMEin molar flow rate of DME at inlet 

FDMEout molar flow rate of DME at outlet 

Fmou molar flow rate of hydrogen at outlet 

Gcs) | molar Gibbs free energy of solid carbon 

Gc partial molar Gibbs free energy of gas carbon 

Gcs) partial molar Gibbs free energy of solid carbon 

AG hus) standard Gibbs function of formation of solid 
carbon 

AG*, standard Gibbs function of formation of species i 

G? standard Gibbs free energy of species i 

Gi partial molar Gibbs free energy of species i 

G' total Gibbs free energy 

ne mole of carbon 

N number of species in the reaction system 

P pressure of system (kPa) 

P? standard-state pressure of 101.3 kPa 

R molar gas constant 

T temperature of system (°C) 

Te temperature at which the first disappearance of 
carbon was achieved (°C) 

I reforming temperature (°C) 

yi gas phase mole fraction 

Greek symbols 

ĝi fugacity coefficient of species i 

Àk Lagrange multiplier 

Wi chemical potential of species i 


of 200-350 °C for MeOH and 200-450°C for DME. DME is 
less toxic and therefore is preferable to MeOH. Reforming of 
DME is expected to match with wide range operating tempera- 
ture of fuel cells. Development of catalysts for DME reforming 
has been currently carried out [6-12]. Durability of the DME 
reforming catalyst is one of the key factors for the practical use. 
Carbon (coke) formation and metal sintering are known as main 
reasons for deactivation. Operating the reforming at appropriate 
condition can suppress such deactivation. Several reports cal- 
culated the thermodynamic composition of reforming processes 
such as MeOH SR [13], methane dry reforming [14], and EtOH 
SR [15] to figure out the condition to avoid the formation of 
coke. Thermodynamic composition of DME SR was recently 
reported [16]. Carbon formation was not however considered in 
the calculation. 

In the present study, thermodynamic equilibrium of 
DME-steam reforming in an external reformer for fuel cell 
applications was studied. Influences of steam-to-carbon ratio 
(S/C = 0-5) and reforming temperature (25—1000 °C) on carbon 


formation boundary, DME conversion, and hydrogen yield were 
evaluated by using Gibbs free minimization technique. 


2. Modeling and simulation methodology 
2.1. Minimization of Gibbs free energy 


Gibbs free energy is the most commonly used function to 
identify the equilibrium state. A minimization of total Gibbs 
free energy is an appropriate method to calculate the equilibrium 
compositions of any reacting system [17]. The total Gibbs free 
energy of a system is given by the sum of ith species: 


N N 
C= SoniGi = Sonim 
i=l i=l 


= X nG? + R(T + 273.15) S— nj Ing (1) 


1 


where G' is the total Gibbs free energy, G; the partial molar Gibbs 
free energy of species i, G? the standard Gibbs free energy, ju; the 
chemical potential, R the molar gas constant, T the temperature 
of system, P the pressure of system, fi the fugacity in system, 
fF the standard-state fugacity, and n; is the mole of species 
i. For reaction equilibria in gas phase, fi = yid;P, i = P®, 
and G? = AG‘, are assumed. By using the Lagrange multiplier 
method, the minimum Gibbs free energy of each gaseous species 
and that of the total system can be expressed as Eqs. (2) and (3), 
respectively: 


:p; P 
AG, + R(T + 273.15) ne +S Arak = 0 (2) 
k 


N = 
ip; P 
Soni (as; + R(T + 273.15) 11216 + re) =0 
i=1 k 
(3) 


with the constraining equation: 
X niar = Ák 
i 


where AGA, is the standard Gibbs function of formation of 


species i, P? the standard-state pressure of 101.3 kPa, y; the gas 
phase mole fraction, 4; the fugacity coefficient of species i, A, 
the Lagrange multiplier, aig the number of atoms of the kth ele- 
ment present in each molecule of species i, and Ax, is the total 
mass of kth element in the feed. 

When solid carbon (graphite) is involved in the system, 
exploiting the vapor—solid phase equilibrium is applied to the 
Gibbs-energy of carbon as shown in Eq. (4). Substituting Eq. (1) 
by Eq. (2) for gaseous species and by Eq. (4) for solid species 
gives the minimization function of Gibbs-energy as following 
Eq. (5): 


Ge) = Ges) = Gcs) = KG Fey =0 (4) 
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yiiP 
k 


HAG g) = 0 (5) 


N-1 


i=1 


ni (as; + R(T + 273.15) In 


where Gog), Gcs), Geis), AG Zog) and ne are the partial molar 
Gibbs free energy of gas carbon, that of solid carbon, the molar 
Gibbs free energy of solid carbon, the standard Gibbs function 
of formation of solid carbon, and mole of carbon, respectively. 

The equilibrium calculations employing the Gibbs-energy 
minimization were done with the Aspen plus, Aspen Tech IM, 
The program is capable to simulate a single phase or multiphase 
of multicomponent in equilibria. The steam-to-carbon ratio S/C 
and the reforming temperature T, were varied in the range of 0-5 
and 25-1000 °C, respectively. The pressure was fixed at 1 atm in 
this study. The equation of state used in the calculation was the 
Peng Robinson method. To perform the calculation, reactant and 
product species with their proportion along with reaction con- 
ditions, i.e. temperature and pressure have to be clarified. Then 
the minimization could be performed to calculate the equilib- 
rium composition. 


2.2. DME SR in external reformer 


The major gas species involved in the DME SR are 
CH30CH3, CH30H, H20, Ho, CO, CO2, CHy and coke, based 
on experimental observations [6-12]. The basis set of com- 
pounds is, therefore, acceptable for practical condition. Steam 
reforming of DME (Eq. (6)) is a two-step reaction, namely, 
hydrolysis of DME to MeOH (Eq. (7)), followed by steam 
reforming of MeOH (Eq. (8)). DME hydrolysis actively takes 
place over acidic sites of acid catalysts such as alumina and zeo- 
lite, while MeOH SR proceeds over metal catalysts, Cu-, Pd-, 
and Pt-based catalysts. Besides DME SR, reverse water gas shift 
reaction (r-WGSR as Eq. (9)) generally proceeds over such metal 
catalysts during the reforming process. In addition, methane can 
be generated via DME decomposition (Eq. (10)) in the case that 
strong acidic catalyst or high reforming temperature is employed 
and via CO methanation (Eq. (11)) in the case Ni-based catalysts 
or precious metal catalysts are used: 


DMESR: (CH3)20 + 3H20 < 6H2+2CO2 (6) 
DMEhydrolysis: (CH3)20 + H20 <> 2CH30H (7) 
MeOHSR: CH30H + H20 <> 3H? + CO2 (8) 
r-WGSR: 0O2+ĦH2< H20 + CO (9) 
DMEdecomposition: (CH3)20 <> CH4 +CO + H2 (10) 
COmethanation: CO + 3H2 <> CH4 + H20 (11) 


Carbon formation boundary was determined by adding car- 
bon (graphite) as a solid form into the basis set. The temperature, 
at which the first disappearance of carbon Te was achieved, was 
considered as a carbon boundary. The regions above and below 
the boundary are the coke-free and coke-formed regions, respec- 
tively. All plausible products DME, MeOH, H20, Ho, CO, CO2, 


CH4, and graphite were taken into account as Case 1. Methane 
formation can be suppressed by using copper catalyst or non- 
acidic catalyst and therefore methane was excluded from the 
calculation as Case 2. Cases 3 and 4 were set for the case when 
carbon formation was not considered. 


e Case 1. CH30CH3, CH30H, H20, Ho, CO, CO2, C 
(graphite), and CH4. 

e Case 2. CH30CH3, CH30H, H20, H2, CO, CO2, and C 
(graphite). 

e Case 3. CH30CH3, CH30H, H20, H2, CO, and CO2. 

e Case 4. CH30CH3, CH30H, H20, H2, CO, CO2, and CH4. 


2.3. DME conversion and hydrogen yield 


To evaluate the performance of the steam reforming system, 
the equilibrium conversions of DME and yield of Ho are defined 
as follows: 


equilibrium conversion of DME (%) 


— FpMEin — FDMEout 


x 100 (12) 
FpMEin 


Ho yield (%) = 


Fetyout [PME 
= 


x 100 13 
ib I, (13) 


FDMEin 


where Fpmgin and Fpmeout are the molar flow rates of DME at 
inlet and outlet, respectively, and Four is molar flow rate of 
hydrogen at outlet. [DME/H2]r = [1/6] is the theoretical mole 
ratio of DME fed and hydrogen produced (see Eq. (6)). 


3. Results and discussion 
3.1. Coke formation boundary 


Formation of coke during the catalytic steam reforming could 
lead to deactivation of catalysts, resulting in low durability and 
activity. Thus it is important to keep it under control. Operating 
the reforming system under the coke-free region could avoid the 
coke formation. We have evaluated the temperature region where 
coke is present or absent in the reforming products. When coke 
with methane or coke is thermodynamically favored, the basis 
sets Casel and Case2 were used for the calculation. Besides CO 
and CO», CH4 was considered as the most abundant carbon- 
containing gaseous product produced during DME SR. 

The most probable reactions leading to carbon formation are 
listed below: 


2CO > CO2+C (14) 
CH4< 2H2+C (15) 
CO +H < HO+C (16) 
CO2 + 2H; > 2H20 + C (17) 
(CH3)20 <> H20 + 2H2 +2C (18) 
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=-= = = = Reaction 14 

Reaction 15 
— = = Reaction 16 
— - —Reaction 17 


— — Reaction 18 
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Fig. 1. Coke formation via reactions (14)—(18) as a function of temperature. 


Fig. | shows the carbon formation as a function of temperature 
via each reaction. This is to determine the role of each reac- 
tion on carbon formation thermodynamically. Reactions (16) 
and (17) are thermodynamically unfavorable at low temper- 
ature than reaction (14) (Boudouard reaction), while reaction 
(15) is favorable at high temperature. Reaction (18) is favorable 
throughout the temperature range studied, since the system is 
not included other carbon-containing compounds, such as, CO, 
CO7, and CH4. 

Fig. 2a and b illustrates the ratio of coke formed per DME 
fed in DME SR as a function of S/C and temperature. In Fig. 2a 
(Case 1), coke was always produced when S/C is less than 0.5. 
Approximately 1 mole of coke was formed by 1 mole of DME 
in the absence of steam. At a given temperature, increase in S/C 
lowered the coke formation. Increase in temperature resulted 
in decrease in coke production except in the temperature range 
of ca. 400-600 “C. However, further increase in temperature 
above 600°C could inhibit coke formation. This is because 
steam reforming and CO, reforming of methane strongly con- 
tributed to the methane removal [18]. In Fig. 2b (Case 2), when 
CH; was exclude from the basis set, coke was always produced 
when S/C is less than 0.5. Two moles of coke was formed by 
1 mole of DME in the temperature range below 200°C, corre- 
sponding to the DME decomposition, reaction (18). Increase 
in temperature and S/C suppressed the coke formation by 
concomitant contribution of reverse reactions (14), (16), and 
(17). 

Coke formation boundary was plotted as a function of S/C 
ratio as shown in Fig. 3. In Case 1, the T, drastically decreased 
from ca. 1100 to 200°C as S/C increased from 0.5 to 1.75. 
Further increase in S/C from 1.75 to 5 gradually decreased Te 
from 200 to 25°C. In Case 2 when methane was not taken 
into account, Te decreased from ca. 900 to 500°C, when the 
S/C ratio increased from 0.5 to 5. The temperature at which 
carbon was unfavorable coincided well with the operating tem- 
perature of MCFC and SOFC. When S/C is lower than 0.5, the 
coke was still present even at high temperature above 2000 °C 
for both cases. When DME SR is operated under the coke- 
formed region, the catalyst that has high resistance to coke 
formation should be used. Considering Case 2 when selectiv- 


Coke formed per DME fed (mol/mol) 


(a) 


Coke formed per DME fed (mol/mol) 


Fig. 2. Coke formation from DME SR as a function of steam-to-carbon (S/C) 
and temperature: (a) Case 1 and (b) Case 2. 


ity to CH4 is zero, the catalyst with high thermal durability is 
required for reforming DME in the coke-free region at high 
temperatures. 


3.2. DME conversion 


Fig. 4 shows the conversion of DME as a function of S/C 
and temperature. In Fig. 4a, DME was completely converted 
throughout the operating condition studied in the presence of 
coke and/or methane (Cases 1, 2, and 4). Thermodynamically 
DME can be converted to methane and coke even at low tem- 
peratures. In Case 2, DME was decomposed via reaction (18). 
When methane is favorable, it could be formed instead of coke 
and hydrogen at the low temperature zone. In Fig. 4b (Case 3), 
when coke and methane are not taken into account, complete 
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Fig. 3. Coke formation boundary of DME SR as a function of steam-to-carbon 
(S/C) and temperature. 


conversion of DME could be achieved when the temperature 
and S/C are higher than 200°C and 1.5 (the stoichiometric ratio 
of DME SR), respectively. Increase in S/C and in temperature 
enhanced the conversion until it reached 100%. However, the 
complete conversion could not be achieved when S/C was below 
0.5. DME could not be reformed in the absence of steam in this 
case. According to the experimental data, DME could be com- 
pletely converted when the reforming temperature was above 
350°C [9-11]. 


3.3. Hydrogen yield 


Fig. 5 shows the hydrogen yield of DME SR as a function of 
S/C and temperature. Theoretically, 1 mole of DME gives 6 mole 
of hydrogen by steam reforming of DME. 

In Fig. 5a (Case 1), H2 yield was promoted by increasing 
temperature. When S/C is higher than 1.5, the maximum yield 
of 90% was found at T, range of ca. 600-700 °C. Further increase 
in T, brought about the decrease in the yield due to the r-WGSR. 
Ho was not produced at the temperature lower than ca. 200°C. 
In Fig. 5b (Case 2), Ho yield was obtained at 33.3% when the 
temperature was below ca. 400 °C. This result is consistent with 
DME conversion and coke formation contributed by reaction 
(18). Hz yield increased when S/C and temperature increased 
up to 5 and ca. 500°C, respectively. The maximum yield of 
95% was achieved at S/C =5 and T, = 500°C. 

As for Case 3 in Fig. 5c, increasing S/C improved the Ho 
yield since DME conversion was expedited and r-WGSR was 
inhibited. The Ho yield of above 99% was achieved at S/C>2 
and T, = 125-250 °C. Upon a complete conversion of DME at a 
given S/C ratio, the Hz yield decreased monotonously as the 
reforming temperature increased. The thermodynamic result 
shows similar tendency with the experimental one although 
the experiment proceeds below the equilibrium [9,11]. High 
space time may be needed for approaching the equilibrium 
of DME SR. As shown in Fig. 5d (Case 4), the Ho yield 
was lower than that obtained from Case 1 at low S/C ratios 
because methane was formed instead of hydrogen and coke. 
When the reforming condition fell into the coke-free region of 


Case 1, 2,4 


DME conversion (%) 
2 


DME conversion (%) 


(b) 


Fig. 4. Equilibrium conversion of DME as a function of steam-to-carbon (S/C) 
and temperature: (a) Cases 1, 2 and 4 and (b) Case 3. 


Case 2, the Ho yield from Cases 1, 2, and 4 becomes compa- 
rable. Even at S/C=0, hydrogen yield could be achieved up 
to 50% yield for Cases 1 and 2, and to ca. 18% for Case 4. 
Hydrogen could not however be produced when the tempera- 
ture were below 200 and 400°C for Cases 1 and 4, respectively. 
Coke formation should be seriously considered under the dry 
condition. 

High S/C and T, enhanced the reforming performance in 
terms of carbon inhibition, DME conversion and hydrogen 
yield. However, the higher S/C requested higher reactor vol- 
ume because of higher steam volumetric flow, and consumed 
higher input heat duty because of higher vaporization energy. 
Therefore, optimization of overall efficiency of fuel cell system 
is required for selection of suitable operating condition of the 
reformer. 
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Fig. 5. Equilibrium hydrogen yield as a function of steam-to-carbon (S/C) and temperature: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4. 


4. Conclusions 


Thermodynamic equilibrium of dimethyl ether steam reform- 
ing (DME SR) was studied by Gibbs free minimization for 
carbon formation boundary, DME conversion and hydrogen 
yield in an external reformer. Effect of steam-to-carbon ratio 
(S/C = 0-5) and reforming temperature (25-1000 °C) and prod- 
uct basis species were investigated. The results indicated that 
carbon formation could be avoided by increasing the steam-to- 
DME ratio and/or by increasing the reforming temperature. 

Based on the compound basis set DME, methanol, CO2, CO, 
Ho, H20 and coke, complete conversion of DME and hydrogen 
yield above 78% were achieved in the coke-free region at the 
normal operating temperature of 600°C for MCFC and that of 
900°C for SOFC. When methane was taken into account, coke 


formation was significantly suppressed. Hydrogen yield up to 
almost 100% could be achieved at S/C > 2 and T, = 125-250 °C 
when coke and methane were thermodynamically unfavorable. 
It should be realized that in the experiment carbon may not be 
formed as graphite and consequently the calculation data may 
deviate from the experiment. Reaction mechanisms, i.e. kinetic 
control or heat and mass transfer control, would also affect the 
result of DME SR in practical applications. 
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